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•  A  model  based  on  dynamic  responses  for  LiFeP04  battery  is  developed. 

•  Effects  of  current  collectors  on  LiFeP04  battery  are  considered  in  this  model. 

•  Dynamic  analysis  for  pulse  behavior  and  electrochemical  rates  are  conducted. 

•  Endothermic  and  exothermic  phenomena  in  different  discharge  rates  are  discussed. 
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An  electrochemical-thermal  model  is  developed  to  predict  electrochemical  and  thermal  behaviors  of 
commercial  LiFeP04  battery  during  a  discharging  process.  A  series  of  temperatures  and  lithium  ion 
concentrations  dependent  parameters  relevant  to  the  reaction  rate  and  Li+  transport  are  employed  in 
this  model.  A  non-negligible  contribution  of  current  collectors  to  the  average  heat  generation  of  the 
battery  is  considered.  Simulation  results  on  rate  capability  and  temperature  performance  show  good 
agreement  with  the  literature  data.  The  behavior  of  Li+  distribution  at  pulse-relaxation  discharge,  the 
variation  of  electrochemical  reaction  rate  and  thermal  behavior  at  a  constant  current  discharge  are 
studied.  Results  of  pulse-relaxation  discharge  describe  the  dynamic  change  of  Li+  concentration  distri¬ 
bution  in  liquid  and  solid  phases,  which  is  helpful  to  analysis  the  polarization  of  the  battery.  In  constant 
current  discharge  processes,  the  electrochemical  reaction  rate  of  positive  electrode  has  a  regular  change 
with  the  time  and  the  position  in  the  electrode.  When  discharge  finished,  there  is  still  a  part  of  the 
LiFeP04  material  has  not  been  adequately  utilized.  At  low  rate,  the  discharge  process  accompanies 
endothermic  and  exothermic  processes.  With  the  rate  increasing,  the  endothermic  process  disappears 
gradually,  and  only  exothermic  process  left  at  high  rate. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Lithium  ion  battery  is  nowadays  one  of  the  most  popular  energy 
storage  devices  due  to  high  energy,  power  density  and  cycle  life 
characteristics  1,2].  It  has  been  known  that  the  overall  perfor¬ 
mance  of  batteries  not  only  depends  on  electrolyte  and  electrode 
materials,  but  also  depends  on  operation  conditions  and  choice  of 
physical  parameters  [3].  Designers  need  an  understanding  on 
thermodynamic  and  kinetic  characteristics  of  batteries  that  is 
costly  and  time-consuming  by  experimental  methods.  Conversely, 
numerical  modeling  and  simulation  for  batteries  are  economic, 
which  can  provide  guidelines  for  design  in  a  short  time  [4],  and 
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information  during  the  electrochemical  and  transport  process.  For 
example,  the  local  electric  potential  and  ion  concentrations 
throughout  the  porous  electrodes  can  be  calculated  while  the  data 
cannot  be  obtained  experimentally  [5,6]. 

The  accuracy  of  numerical  modeling  and  simulation  of  electro¬ 
chemical  and  thermal  behavior  relies  on  the  model  construction 
and  the  parameters  applied  during  simulation  [7].  The  most  famous 
and  practical  model  for  lithium  ion  battery  is  the  porous  electrode 
model  [8,9],  which  was  based  on  the  porous  electrode  theory 
containing  charge  transfer  kinetics  at  reaction  sites,  species  and 
charge  conservations.  It  was  combined  with  an  energy  conserva¬ 
tion  equation  by  Newman  and  Pals  [  10,11  ,  and  was  developed  as  an 
electro-thermal  model,  which  made  it  possibly  simulate  and  pre¬ 
dict  the  interaction  between  temperature  and  electrochemical  re¬ 
action  of  batteries.  Previously,  the  electrochemical  models  and 
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Nomenclature 

t+ 

Li+  transference  number 

T 

absolute  temperature  (I<) 

Lamb 

ambient  temperature  (I<) 

List  of  symbols 

Ui 

open  circuit  voltage  (V) 

Acell 

area  of  the  positive  electrode  (both  sides)  (m2) 

V 

thermodynamic  factor  relating  to  electrolyte  activity 

Cl,i 

lithium  in  active  material  (mol  m“3) 

X 

distance  from  half  negative  foil  along  negative-positive 

Cl, max, i 

maximum  concentration  (mol  m-3) 

direction  (m) 

^1, surf, i 

Li+  concentration  on  the  surface  of  active  material 

particles  (mol  nrr3) 

Greek  letters 

CP,i 

heat  capacity  (J  (kg  K)_1) 

«  a,i 

transfer  coefficient  for  anodic  current 

Du 

solid  phase  diffusivity  (m2  s-1) 

0LCJ 

transfer  coefficient  for  cathodic  current 

£ho,i 

reference  solid  phase  diffusivity  (m2  s-1) 

£l,i 

active  material  volume  fraction 

£a,D,i 

diffusion  active  energy  (kj  mol-1) 

£2,i 

volume  fraction 

£a,k,i 

reaction  active  energy  (kj  mol-1) 

(Pi 

electric  potential  (V) 

h 

heat  transfer  coefficient  (W  nr2  K-1) 

Ti 

Bruggeman  exponent 

I 

current  (A) 

k 

ionic  or  electronic  conductivity  (S  nr1) 

Japp 

cell  current  density  related  to  Aceii  (A  m-2) 

Pi 

density  (kg  nr3) 

Jo.i 

exchange  current  density  (A  m-2) 

(Li 

solid  phase  conductivity  (S  nrr1) 

Jloc.i 

local  current  density  (A  m-2) 

fco.i 

reaction  rate  constant  (m25  mol-0-5  s_1) 

Subscripts  and  superscripts 

ki 

thermal  conductivity  (W  (m  K)-1) 

0 

initial  or  equilibrated  state 

U 

thickness  of  component  (m) 

1 

solid  phase 

Qact 

active  heat  generation  (J  m-3) 

2 

liquid  phase 

Qohm 

ohmic  heat  generation  (J  m-3) 

amb 

ambient  temperature 

Qrea 

reaction  heat  generation  (J  m-3) 

n 

negative  electrode 

r 

radius  distance  variable  of  particle  (m) 

P 

positive  electrode 

Ri 

characteristic  radius  of  electrode  particles  (m) 

irr 

irreversible 

5a,  i 

specific  surface  area  (nr1) 

re 

reversible 

SOC0,i 

initial  state  of  charge 

s 

separator 

t 

time  (s) 

electrochemical  thermal  models  were  built  for  simulating  the 
cathode  materials  of  lithium  ion  battery,  such  as  LiCo02,  LiMn204  or 
Li(NiCoMn)02  [12-18].  Later,  based  on  the  observation  of  a  phase 
change  occurred  in  the  LiFeP04  cathode  during  the  lithiated  and 
unlithiated  process  [19,20],  Srinivasan  [21]  developed  a  model  that 
accounted  for  the  phase  change  with  the  shrinking  core,  and 
investigated  the  cause  for  the  low  power  capability  of  the  materials. 
The  other  models  without  special  features  about  the  two-phase 
process  were  also  founded  for  further  qualitative  analysis  [22- 
31].  For  instance,  Wang  [22]  analyzed  the  effect  of  local  current 
density  on  electrode  design  for  the  LiFeP04  battery.  Ye  [23]  devel¬ 
oped  an  electro-thermal  cycle  life  model  by  incorporating  the 
dominant  capacity  fading  mechanism  to  investigate  the  capacity 
fading  effect  on  the  performance.  Gerver  [24]  and  Christian  Flellwig 
[25]  gave  a  multidimensional  modeling  framework  for  simulating 
coupled  thermal  and  electrochemical  phenomena  that  are  critical 
for  safety,  durability  and  design  optimization  studies.  These  models 
demonstrated  that  the  porous  electrode  theory  can  be  used  to 
simulate  a  LiFeP04  battery  system  without  considering  the  phase 
change. 

In  the  previous  work  [8-31],  the  computational  domains  were 
divided  into  three  parts,  named  negative  electrode,  separator  and 
positive  electrode,  but  the  current  collectors  of  a  battery  were 
neglected.  Because  a  rate  of  the  heat  generation  is  an  average  value 
obtained  through  dividing  the  total  heat  generation  by  the  cell 
volume,  about  10%  of  the  cell  volume  comes  from  the  current  col¬ 
lectors,  indicating  that  the  contribution  of  the  collectors  to  the  heat 
generation  rate  cannot  be  neglected. 

Besides,  the  material  properties  treated  as  parameters  also  have 
an  important  influence  on  the  accuracy  of  the  simulation  [7]. 
During  the  discharge  process,  lithium  ions  deintercalate  from  the 
negative  porous  electrode,  then  transfer  through  separator  and 


intercalate  into  the  positive  electrode,  resulting  in  changes  of 
temperature  and  lithium  ion  concentrations  in  the  solid  and  liquid 
phases  that  may  affect  the  diffusion  coefficient,  ionic  conductivity, 
ion  transference  number  as  well  as  the  reaction  rate  constant 
[24,32-34].  Therefore,  it  is  necessary  to  rectify  dynamically  these 
parameters  during  a  simulation. 

In  this  study,  we  develop  an  electrochemical  thermal  model 
for  a  LiFeP04  battery  by  considering  the  current  collectors  into 
the  computational  domain,  dynamic  responses  in  lithium  ion 
concentration,  and  temperature  as  parameters  during  the 
discharge  process.  By  comparing  experimental  results  with 
simulation  at  different  operating  temperatures  and  discharge 
rates,  this  model  can  be  used  to  study  the  dynamic  evolution  for 
pulses,  relaxation  behavior,  electrochemical  reaction  and  thermal 
behavior  at  a  constant  discharge  rate  in  lithium  iron  phosphate 
battery. 

2.  Model  development 

2.1.  Model  assumption  and  simulation  domain 

This  electrochemical— thermal  model  for  a  LiFeP04  battery  is 
developed  based  on  the  porous  electrode  model  [8,9].  The  active 
materials  of  solid  electrodes  are  treated  as  homogenous  media,  and 
are  comprised  with  spherical  particles.  Since  the  current  collectors 
play  an  important  role  in  electron  and  heat  transferring,  the  current 
collectors  are  considered  as  a  computing  domain. 

Fig.  1  shows  a  schematic  computational  domain  of  one¬ 
dimensional  (ID)  battery  model  [8,9].  Because  both  positive  and 
negative  current  collectors  have  two  sides  contacting  and  reacting 
with  electrolyte,  we  take  a  half  of  their  thickness  into  simulation 
domain  in  order  to  accurately  calculate  the  current  flow  and  the 
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Fig.  1.  Schematic  diagram  of  lithium  iron  phosphate  battery  and  computational  domain. 


heat  generation.  There  are  six  boundaries  in  the  model.  Four  inner 
boundaries  (anode  current  collector/anode  interface  boundaries  2, 
anode/separator  interface  boundaries  3,  separator/cathode  inter¬ 
face  boundaries  4  and  cathode/cathode  current  collector  interface 
boundaries  5)  and  two  external  boundaries  (anode  current  collec¬ 
tor  external-face  boundary  1  and  cathode  current  collector 
external-face  boundary  6)  are  shown  in  Fig.  1. 

2.2.  Electrochemical  part 

2.2.1.  Electronic  charge  conservation 


14s 


2RT 

-  V(/>2  +  “^" 


1  + 


9ln  / 


din  c2_ 


(1-M^ 

c2  . 


kf  =  k2e\- 


*^a,i/loc,i 


(6) 

(7) 


Liquid-junction  potential  is  introduced  in  Eq.  (6)  with 
expression 


^junc  — 


2  RT 

~r 


1  + 


ain  / 


Oln  c2_ 


M  ,  \  2 RT 

a -4)  =  -p-v 


(8) 


(a)  Solid  phase.  There  is  a  double  layer  capacitance  at  the 
interface  between  the  active  materials  and  electrolyte  [35  .  The 
capacitance  can  store  additional  energy  and  smooth  the  abrupt 
change  of  the  cell  voltage  caused  by  the  short-time  pulses  passing 
through  a  cell.  However,  the  capacitance  was  seldom  considered  in 
the  previous  models  [35].  Small  modification  in  the  electronic 
charge  conservation  is  required  to  account  of  the  double  layer. 
Electronic  charge  conservation  for  solid  phase  can  be  expressed  as 
follows: 


The  parameter  v  is  the  thermodynamic  factor  relating  to  elec¬ 
trolyte  activity,  and  it  is  concentration  and  temperature  dependent. 

There  is  no  flux  at  external  boundaries  (boundary  1  and 
boundary  6),  which  is  expressed  as  Eq.  (9);  the  parameter  cp2  is 
taken  to  be  continuous  at  boundary  2  and  boundary  5. 


0X 


x=L 


ncc 


0X 


X—Lncc  +tn  +tS  +tp 


=  0 


(9) 


2.2.2.  Mass  conservation 


V-(-(7cV0c)  =  -Ji  (1) 


where  cpc  is  arbitrarily  set  to  zero  at  boundary  1 ;  at  boundary  6,  the 
charge  flux  is  set  to  equal  the  average  current  density  of  the  battery, 
which  is  expressed  as  Eq.  (4);  at  boundary  3  and  boundary  4,  there 
is  no  charge  flux,  which  is  expressed  as  Eq.  (5)  and  the  boundary 
condition  is  set  as  isolation. 

^ c lx=0  =  _<7cV0c|x=incc+Ln+Ls+Lp+Lpcc  =  -/app  (4) 

-kfvfa  I  =  0; I  =0  (5) 

'X—. Lncc+Ln  lx— Lncc+Ln+Ls 


(b)  Solution  phase.  According  to  the  concentrated  solution 
theory,  the  governing  equation  for  electronic  charge  conservation 
in  solution  phase  is  expressed  as: 


(a)  Solid  phase.  The  mass  conservation  of  lithium  ions  in  an 
intercalation  particle  of  the  electrode  active  material  is  described 
by  Fields  law.  The  parameter  r  is  the  distance  from  the  center  of 
solid  sphere,  the  mass  transport  within  solid  phase  in  spherical 
coordinates  can  be  described  as: 


9ci,i 

dt 


I"9  ! 

^r20Cl,A 

dr ' 

V  dr  )\ 

=  D 


l,i 


2dcu 
r  dr 


(10) 


Let  the  variable  t  equal  to  r/R[ ,  the  Eq.  (10)  can  be  described  as 


T2p9Cl,i  Dl.i 

T  Ri~dF  =  R, 


I"9  1 

^T20Cl,iY 

9t  ' 

V  dz  )\ 

(11) 


The  parameter^  is  the  particle  radius.  There  is  no  species  source 
at  the  center  of  sphere,  so  0Citi/0T|T=o  =  0.  The  Li+  concentration  on 
the  surface  of  the  particles  is  coupled  to  the  concentration  and  flux 
in  the  ID  model  for  the  charge  and  material  transport  in  the 
electrolyte. 


(b)  Solution  phase.  Solution  phase  mass  conservation  for  LiF6 
dissolved  in  the  liquid  phase 
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^  +  V-{-DfVC2}  ^(1-4) 
Df  =  D2el> 


(12) 

(13) 


Reaction  heat  generation  is: 


Qrea 


dU{ 

dr 


Ohmic  heat  generation  is: 


(24) 


The  flux  of  liquid  species  is  arbitrarily  set  to  zero  at  boundary  2 
and  boundary  5,  liquid  species  flux  and  species  concentration  at 
boundary  3  and  boundary  4  are  taken  to  be  continuous. 


2.2.3.  Electrochemical  kinetics 

The  local  current  per  active  material  area  is  calculated  using  the 
Butler-Volmer  equation: 


Jloc,i  -  J 0,ii  exP 


«a,i  ViF' 


RT 


-  exp 


-«C,  iVif 


RT 


(14) 


The  parameter  jioc,i  is  driven  by  over-potential,  the  parameter  rj[ 
is  defined  as  the  difference  between  solid  and  electrolyte  phase 
potentials  minus  Ui. 

Vi  =  01, i  -  02, i  -  Ui  (15) 


The  parameter  Ui  is  the  thermodynamic  equilibrium  potential  of 
the  solid  phase  and  is  taken  to  be  a  function  of  the  local  SOC  on  the 
surface  of  active  particles. 

The  temperature-dependent  open  circuit  potentials  of  positive 
and  negative  electrodes  are  approximated  by  Taylor’s  first  order 
expansion  around  a  reference  temperature: 

Uj  =  (Jlef,i  +  (t  -  Tref )  ^  (16) 

The  parameter  L/ref,i  is  the  open  circuit  potential  under  the 
reference  temperature. 

In  Eq.  (14),  exchange  current  density,  jo, i,  connects  concentra¬ 
tions  in  both  solid  and  liquid  phase: 

Jo,i  =  Fk[c2  (cl,max,i  —  cl,surf,i)  Cl,surf,i  (^) 


The  parameter  k\  is  the  reaction  rate,  considered  temperature 
dependent  in  this  paper.  The  parameters  aa  and  ac  are  the  anodic 
and  cathodic  transfer  coefficients,  respectively. 

The  outputs  of  the  model  are  the  cell  potential,  current  density 
distribution,  species  and  concentrations  distributions,  and  the  cell 
potential  is  derived  by  following  expression: 

£  =  01  lx=LnCc+In+Is+Ip+Ipcc  ~  ^1  lx=0  H8) 

2.3.  Energy  conservation 

The  heat  generation  is  modeled  with  the  local  heat  generation 
model  of  Rao  and  Newman  [36]  using  the  formulation  of  Gu  and 
Wang  [37].  According  to  the  position  of  heat  generation,  the  total 
heat  generation  is  the  summation  of  heat  generated  in  the  two 
electrodes,  separator  and  current  collectors.  The  summary  of  the 
heat  generation  mechanisms  is  listed  in  Table  1. 

According  to  the  type  of  the  heat  generation,  there  are  three 
parts  of  heat  sources  during  the  charge  or  discharge  processes, 
including  reaction  heat,  Qrea,  due  to  entropy  change  during  the 
discharge;  ohmic  heat,  Qohm,  due  to  the  ohmic  potential  drop  and 
active  polarization  heat,  Qact,  due  to  the  electrochemical  reaction 
polarization  between  active  material  particle  surface  and  the 
electrolyte. 


Qohm  —  Qi,(2)  +  Qi,(3)  +  Qs  +  Qi,c  (25) 

Active  heat  generation  is: 

Qact  —  Sa,ijloc,i(01,i  —  02,  i  —  ^i)  (26) 

Active  heat  Qact  and  ohmic  Qohm  are  irreversible,  while  reaction 
heat  Qrea  is  reversible. 

Qirr  =  Qohm  +  Qact  (27) 

Ore  —  Qrea  (28) 

The  energy  conservation  in  lithium  ion  battery  is  shown  as 
follows: 

ftCp,i^  +  V-(-kiVT)  =  Qjrr  +  Qre  (29) 


The  parameters  pi,  Cp,i  and  /q  are  density,  heat  capacity  and 
thermal  conductivity,  respectively. 

According  to  Newton’s  cooling  law  and  radiation  law,  the 
boundary  condition  for  energy  conservation  is  expressed  as: 

-Avr  =  -fi(Tamb  —  T)  —  eo  (ramb  -  r4)  (30) 

The  parameter  A  is  thermal  conductivity  of  stainless  steel  case,  h 
is  natural  convection  heat  transfer  coefficient  (h  =  7.17  W  m-2  K-1 
[38]),  Tamb  is  ambient  temperature,  and  e  is  the  blackness  of  the 
battery  surface  (e  =  0.8  [38]). 

3.  Model  parameters  and  model  validation 

3.1.  Battery  parameters  and  thermal  properties 

The  physical  properties  of  battery  components  and  battery 
design  parameters  are  summarized  in  Tables  2  and  3,  respectively. 

3.2.  Dynamic  parameters 

The  dynamic  response  of  battery  properties  is  attributed  to  the 
change  of  temperature  and  lithium-ion  concentration  during  the 
discharge  and  charge  processes.  It  is  considered  for  eight  sets  of 
physical  properties.  These  are  (1 )  /<o,i,  the  reaction  rate  constant;  (2) 
Dij,  the  diffusion  coefficient  of  Li+  in  the  solid  active  particles;  (3) 
open  circuit  potential  (OCP)  of  the  electrodes;  (4)  D2,  the  diffusion 
coefficient  of  Li+  in  the  electrolyte;  (5)  /<2,  the  electrolyte  ionic 
conductivity;  (6)  v,  the  thermodynamic  factor  relating  to  electrolyte 
activity;  (7)  t+,  the  Li+  transference  number;  (8)  ac,  current  col¬ 
lector  electrical  conductivity. 

3.2 A.  Dynamic  variables  related  to  electrodes 
3.2  A  A.  Electrode  kinetics  properties.  The  dynamic  temperature  de¬ 
pendences  of  Diti  follow  the  Arrhenius  equation.  However,  the  LFP 
electrode  exhibited  significantly  higher  utilization  on  charging  than 
on  discharging  at  the  same  current  density  [39].  This  asymmetry 
was  ascribed  to  the  different  transport  limitation  between  the  fully 
charged  and  discharged  states  [40].  To  account  for  this  asymmetry, 
an  empirical  equation  is  used  for  the  solid-state  Li+  diffusion  co¬ 
efficient  that  is  dependent  on  the  average  lithium  ion  concentration 
of  particles  [27]: 
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Table  1 

Heat  generation  mechanism  in  a  cell. 
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Heat  generation  mechanism 

Porous  electrode 

(1).  Electrochemical  interface 


(2) .  Electrical  ohmic  heat 

(3) .  Ionic  ohmic  heat 


Separator 
Ionic  ohmic  heat 


Current  collector 
Electrical  ohmic  heat 


Equation 


9i,(l)  _  ^a,i/loc,i  ^  3"  01,i  02,i  t/j ^ 

<?i,(2)  =  4-V014 

2  RTke^ 

Vi, (3)  =  /<f  V02;i-V02)i  +  — p^(4  -  l) 


(19) 


(20) 


1  + 


9ln/ 


aln  c2;i 


•V(ln  c2  j)  •  V<^2i  (21) 


2 RTkeff  r  aln  f  1 

qs  =  Jcfv</>2-V02+  — ^-(4-1)  1+gjj^  -V(lnc2)-V02 


(22) 


Qi,c  —  ^i,c^0i,c '  ^0i,c 


(23) 


Subscript  i  =  n  or  p  for  a  negative  or  a  positive  electrode. 


Dip 


1.18  x  IQ"18 

P+y)1'6 


(31) 


In  this  study,  we  take  a  combination  of  the  Arrhenius  equation 
and  Eq.  (31)  into  the  formula  for  the  D i)P,  which  is  expressed  as: 


D 


i,p  = 


1.18  x  10-18 

P+y)1-6 


exp 


£yp.p  P  1  + 

R  \T  298.15 JJ 


(32) 


The  state  of  charge  (SOC)  and  temperature  also  affect  /<o,i  and  D itn, 
but  duo  to  lack  of  data,  the  SOC  is  not  considered.  I<n  and  Di>n  are 
expressed  by  Eqs.  (33)-(35)  followed  the  Arrhenius  formula  [41  ]. 


The  parameters  Ei, d)P  and  EitD,n  are  the  activation  energy  for 
diffusion  of  lithium  respectively  in  the  positive  and  negative  elec- 
trades,  (Ei,d,p  =  35  kj  mol-1  [23,28],  Ei,D,n  =  35  kj  mol-1  [23,42]). 
The  parameters  E y<tP  and  E\^n  are  the  activation  energy  for  the  rate 
constant  of  the  positive  and  negative  electrodes,  respectively 
(Ei,k,p  =  20  kj  mol-1  [23,43],  E,,k,n  =  30  kj  mol”1  [23,28]). 


3.2A.2.  Electrode  thermodynamics  properties..  The  open  circuit  po¬ 
tentials  of  the  positive  electrode  Up  and  the  negative  electrode  Un 
are  a  function  of  local  SOC  on  the  surface  of  active  particles.  Up  and 
Un  come  from  Ref.  [28],  which  are  expressed  by  Eqs.  (36)  and  (37). 


Di,n  =  3.9xl0-14exp(  — __!_))  (33) 

4*  = 3  x  10 11  exp(-nr(T-2 90s))  (34) 

k^p  =  l.4xl0-12exp(-3y)exp(-%P(l-29TT5)) 

(35) 


Table  2 

Thermal  properties  of  battery  components. 


Materials 

Density 

(kg  ITT3) 

Heat 
capacity 
a  (kg  K)-1; 

Thermal 

conductivity, 

)  (W(mK)-1) 

Electrical 

conductivity, 

(S  ITT1) 

Negative 

2223 

641 

1.04 

100 

electrode 

Positive 

1500 

800 

1.48 

0.5 

electrode 

Separator 

900 

1883 

0.5 

- 

Electrolyte 

1210 

1518 

0.099 

Eqs.  (42)  and  (43) 

Copper  foil 

8700 

396 

398 

-0.04889 T3  +  54.65T2 
-  218.00T  +  3.52  x  106 
(s  cm-1) 

Aluminum 

foil 

2700 

897 

237 

-0.0325 T3  +  37.07 T2 
-  15,0007+2.408 
x  106  (s  cm-1) 

Stainless 

7500 

460 

14 

- 

steel  case 

l/ref  =3.4323  -  0.4828  exp(  -  80.2493(1  -y)1  3198) 

-  3.2474  x  10“6  exp (20.2645(1  -y)3'8003)  (36) 

+  3.2482  x  10~6  exp (20.2646(1  -y)3  7995) 


l/ref  =0.6379 +  0.541 6  exp(-305.5309x) 


+  0.044  tanh 


/  x- 0.1958\ 
V  0.1088  ) 


-0.1978  tanh 


x  —  1.0571\ 
0.0854  ) 


-0.6875  tanh 


x  —  0.01 17\ 
0.0529  ) 


-0.0175  tanh 


x-0.5692\ 
0.0875  ) 
(37) 


Entropy  change  in  electrodes  is  AS  =  nF(dU/dT).  dUp/dT  and 
dUnldl  are  the  entropy  changes  of  lithium  iron  phosphate  positive 
electrode  and  the  negative  electrode,  respectively.  Their  curves 
[44,45]  are  as  shown  in  Fig.  2,  and  are  expressed  by  Eqs.  (38)  and 
(39). 


=  -0.35376y8  +  1.3902y7  -  2.2585y6  +  1.9635y5 

-  0.98716y4  +  0.28857y3  -  0.046272y2  +  0.0032158y 

-  1.9186  x  1CT5 


(38) 
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Critical  battery  parameters  used  in  baseline  for  a  2.3  Ah  LiFeP04  cylindrical  26650  type  battery. 
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Quantity 


Negative  electrode  Positive  electrode  Separator  Electrolyte  Copper  foil  Aluminum  foil 


Design  specifications  (geometry  and  volume  fractions) 


^cell(m2) 

0.1694 

£l,i 

0.55 

0.43 

- 

£2,i 

0.33 

0.332 

0.54 

Li  (|im) 

34 

70 

25 

6.2  (half  thickness) 

10  (half  thickness) 

Ri(\im) 

0.0365 

3.5 

Lithium  ion  concentrations 

Cini,2  (mol  ITT3) 

1200 

Cmax  (mol  m  3) 

31,370 

22,806 

Cini.l  (mol  ITT3) 

31,370  *  0.86 

22,806  *  0.022 

Kinetic  and  transport  properties 

Oa,i>  ^c.i 

0.5 

0.5 

Yi 

1.5 

1.5 

1.5 

D2  (m2  s"1) 

Eq.  (40) 

Cdi  (F  m"2) 

0.2 

0.2 

Di  i  (m2  s"1) 

Eq.  (33) 

Eq.  (32) 

k\  (m2-5  mol-0  5  s"1) 

Eq.  (34) 

Eq.  (35) 

Fa,k,i  (J  mol'1) 

20,000 

30,000 

kl  (S  ITT1) 

100 

0.5 

k2  (S  m_1) 

Eqs.  (42)  and  (43) 

ac 

t+ 

Eq.  (44) 

V 

Eq.  (41) 

Constant  quantity 

rref(T) 

298.15 

F  (C  mol1) 

96,487 

^  =  344.1347148 

dr 

exp(-32.9633287x  + 8.316711484) 

X  1+749.0756003  exp(-34.79099646x  + 8.887143624) 

-  0.8520278805*  +  0.362299229X2  +  0.2698001697 

(39) 

3.2.2.  Dynamic  variables  related  to  lithium  ion  transport  in 
electrolyte 

Due  to  lack  of  the  proportion  of  the  electrolyte  components,  the 
dynamic  variables  dependent  on  temperature  of  lithium  ion 
transport  in  LiPF6  are  assumed  to  follow  the  properties  of  similar 
electrolyte  systems.  Valoen  [32]  investigated  a  full  set  of  transport 
properties  for  LiPF6  in  PC/EC/DMC  (10:27:63  by  volume)  experi¬ 
mentally,  showing  the  k 2,  D2,  and  v  as  functions  of  temperature  and 
LiPF6  concentration. 

D2  =  1  X  10“4  X  1 0  4-43  7~— 229.0— dose  2.2x io~4C  (40) 


v  =0.601  -0.24Vl0-3C 

/  , - N  (41) 

+  0.982(3  -  0.0052(T  -  294.0)  V10-9C3) 

k2  =  lx  10~4(-  10.5  +  0.0747-6.69  x  10-572  +  6.68 
x  10“4C-  1.78  x  10-5CT  +  2.8  x  10-8CT2 
+  4.49  x  10~7C2  -  8.86  x  10-10C2t) 

(42) 

By  combined  the  parameters  of  LiPF6  in  EC/DMC  (2:1  by  volume) 
at  25  °C  [46]  with  that  in  the  PC/EC/DMC  system  [32],  the  following 
expression  is  obtained: 


k2  =  1.12  x  10_4(  — 8.2488  +  0.0532487  —  2.9871  x  10“572 
+  0.26235C- 9.3063  x  10-3C7  +  8.069 
x  10-6CT2  +0.22002C2  -  1.765  x  10“4C27) 

(43) 

The  Li+  transference  number,  t+,  is  considered  as  a  constant, 
such  as  t+  =  0.363  [14].  Flowever,  t+  also  appears  to  be  dependent  of 
Li+  concentration  and  temperature  [47,48].  Hence,  t+  is  adopted 
from  the  published  literature  [47,48]  that  is  a  function  of  temper¬ 
ature  and  Li+  concentration  in  this  study. 

t*,2.67xl0-exp(?|5)(45)+3.09 

10  5  e+6f )(,«»,)  0.517«p(  ®6]  (44) 

It  should  be  noted  that  the  solvent  mixture  for  the  commercial 
battery  is  unknown.  We  assume  that  Eqs.  (40)— (44)  are  applicable 
in  commercial  batteries,  so,  these  data  are  used  in  this  simulation. 

3.3.  Model  validation 

The  electrochemical— thermal  model  is  validated  from  electro¬ 
chemical  performance.  Experiment  data  are  from  the  introduction 
document  of  ANR26650ml-a  lithium  ion  cylindrical  cell  (A123 
company,  US)  [49]. 

Firstly,  the  experiment  data  [49]  of  battery  potential  against  the 
discharge  capacity  are  compared  with  simulation  results  that  are  in 
fairly  good  agreement  with  the  experimental  results  (see  Fig.  3).  In 
order  to  further  validate  the  electrochemical  performance  at 
different  discharge  rates,  the  simulation  results  are  compared  with 
the  experiment  data  [25]  at  different  discharge  rates  (C/10, 1/2C,  1C, 
2C)  at  25  °C  in  Fig.  4.  The  results  indicate  this  electrochemical- 
thermal  model  based  dynamic  response  is  reliable  to  simulate  the 
discharge  performance  of  lithium  iron  phosphate  battery  at 
different  discharge  rates. 


136 


J.  Li  et  al.  /  Journal  of  Power  Sources  255  (2014)  130-143 


Fig.  2.  The  entropy  change  of  the  positive  electrode  and  negative  electrode  as  a 
function  of  SOC. 


4.  Results  and  discussion 

This  model  is  mainly  used  to  predict  the  electrochemical  per¬ 
formance  and  thermal  behavior  of  lithium  ion  batteries. 

4.1.  Dynamic  evolution  for  pulse  behavior 

Lithium  ion  batteries  have  been  used  in  many  aspects  which 
work  intermittently,  for  example,  assisting  the  engine  of  hybrid 
electric  vehicles  during  vehicle  acceleration.  Because  these  actions 
are  of  short  duration,  on  the  order  of  seconds,  the  pulse  and 
relaxation  behaviors  are  considered  important  aspects  to  analyze 
polarization  in  battery.  In  this  study,  pulse  test  with  different 
discharge  rates  was  chosen  to  obtain  Li+  concentration  distribution. 
During  simulated  operation,  the  battery  is  discharged  at  0.2C,  1C,  or 
2C  constant  current  during  0-200  s  and  is  relaxed  by  cutting  off 
current  during  200-1000  s.  The  Boundaries  2, 3, 4  and  5  are  used  to 
examine  the  spatial  and  temporal  distribution  of  lithium  ion  con¬ 
centration  in  liquid  and  solid  phases. 

The  spatial  and  temporal  distribution  of  lithium  ion  concen¬ 
tration  in  electrolyte  (C2)  during  the  pulse  test  at  different  discharge 
rates  is  depicted  in  Fig.  5.  At  the  beginning  of  the  discharge,  a 
distribution  of  lithium  ion  rapidly  builds  up  inside  the  battery,  C2 
decreases  from  Boundary  2  to  Boundary  5,  forming  a  gradient 
which  drives  lithium  ions  along  the  battery  direction  from  the 
negative  electrode  to  positive  electrode.  At  a  high  discharge  rate  (1C 
or  2C),  C2  increases  in  the  negative  electrode  and  decreases  in  the 


positive  electrode  during  the  very  early  period  (about  80  s),  and 
then  stays  steady  from  80  s  to  200  s.  On  the  contrary,  at  a  low 
discharge  rate,  C2  changes  from  0  s  to  200  s.  When  the  current  is  cut 
down  at  200  s,  C2  changes  rapidly  in  a  timescale  of  millisecond. 
Then  a  self-adjusting  period  driven  by  diffusion  is  taken  to  attain 
equilibrium.  The  greater  discharge  rate  is,  the  longer  self-adjusting 
period.  Due  to  the  employment  of  double  layer  capacitance  of  the 
active  materials  and  electrolyte  interface,  C2  in  equilibrium  shows 
an  increase  compared  to  the  initial  level  (1200  mol  m-3).  At  0.2, 1 
and  2C  discharge  rates,  the  augmentation  is  0.6,  8.4,  9.9  mol  m-3, 
respectively,  showing  a  rising  trend  along  with  an  increase  of  the 
discharge  rates.  However,  this  rising  is  not  obvious  at  high 
discharge  rates. 

Lithium  ion  concentration  gradient  in  solid  phase  can  be 
described  as  an  indicator  of  polarization  in  a  particle.  Large  con¬ 
centration  gradient  in  a  solid  particle  means  high  polarization, 
which  may  lead  to  high  diffusion  induced  stress  (DIS),  cause  par¬ 
ticle  fracture  and  reduce  battery  capacity  and  battery  power  [50]. 
As  introduced  in  the  section  of  model  development,  active  mate¬ 
rials  (LiFePCH  and  Graphite)  in  solid  electrodes  are  considered  to  be 
homogenous,  and  they  are  composed  of  spherical  particles.  In  order 
to  study  their  interior  lithium  ion  concentration  gradient,  four 
special  particles  at  Boundary  2, 3, 4  and  5,  respectively,  are  selected, 
and  they  divided  into  10  parts  along  the  radial  direction.  The  center 
and  surface  of  the  particles  are  located  at  r/R  =  0  and  1,  respectively. 
Fig.  6  gives  lithium  ion  concentration  distribution  throughout  the 
particles. 

At  all  of  the  boundaries,  lithium  ion  concentrations  of  the  par¬ 
ticles  show  great  changes  at  pulse  time  of  0-200  s.  They  decrease  at 
Boundary  2  and  3  in  the  negative  electrode  with  discharge  time,  as 
shown  in  Fig.  6a  and  b.  However,  an  increase  with  discharging  time 
at  Boundary  4  and  5  in  the  positive  electrode  is  seen  in  Fig.  6c  and  d. 
At  200  s,  the  surface  concentrations  of  particles  at  Boundary  2  and  3 
decrease  to  25.2  kmol  nrr3  and  24.3  kmol  m-3,  from  initial  con¬ 
centration  (26.98  kmol  m-3),  and  at  Boundary  4  and  5  the  con¬ 
centrations  increase  to  2.3  kmol  m-3  and  1.7  kmol  m-3  from  initial 
concentration  (0.5  kmol  nrr3).  These  indicate  the  lithium  ion  dis¬ 
tribution  is  in  location  dependent. 

After  the  current  is  cut  down  at  200  s,  the  battery  is  in  relaxa¬ 
tion.  In  the  relaxation  period  (200-1000  s),  the  variation  of  lithium 
ion  concentration  in  the  negative  electrode  including  Boundary  2 
and  3  appears  the  similar  tendency.  At  r/Rn  =  0, 0.1, 0.2,  0.3,  0.4,  the 
concentrations  decrease  gradually  until  unchanged.  And  at  r/ 
Rn  =  0.5,  0.6,  0.7,  0.8,  0.9,  1.0  the  concentrations  turn  to  increase 
until  unchanged.  As  for  the  positive  electrode,  the  variation  of 


Capacity  (Ah) 


Fig.  3.  -20  °C,  0  °C,  25  °C,  45  °C,  1C  discharge  validations. 
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interior  lithium  ion  concentration  of  the  particle  at  Boundary  5 
shows  a  contrary  location  dependent  with  Boundary  2  and  3.  At  r/ 
ftp  =  0,  0.1,  0.3,  0.2,  0.4,  the  concentrations  continue  increasing 
gradually  until  unchanged,  and  at  r/ftp  =  0.5, 0.6, 0.7, 0.8, 0.9, 1.0  the 
concentrations  turn  to  decrease  until  unchanged.  But  for  Boundary 
4,  lithium  ion  concentrations  of  the  particle  decrease  gradually 
until  unchanged,  appearing  no  location  dependence  that  is  quite 
different  from  at  other  boundaries. 

In  the  relaxation  period,  the  interior  concentrations  remain 
unchanged  ultimately,  which  is  in  equilibrium.  It  is  well  known 
that  the  cells  of  a  battery  require  a  certain  time  to  reach  a  steady 
state  in  terms  of  charge,  concentration  and  temperature.  As 
introduced  in  literature  [51],  relaxation  time  factor  (ft2/D)  denotes 
the  theory  time  needed  for  interior  concentrations  to  obtain 
equilibrium  after  pulse  discharge,  which  reflects  the  depolariza¬ 
tion  ability  of  the  battery.  In  theory,  at  25  °C,  ftn/Dn  =  136  s,  and 


ftp/Dp  =  1129  s.  With  the  discharge  proceeding,  battery  tempera¬ 
ture  will  increase  that  leads  to  D  increase  (see  Eqs.  (32)  and  (33)). 
So  ft2/D  at  discharge  process  is  lower  than  that  at  25  °C.  That 
explains  why  the  relaxation  times  of  negative  and  positive  elec¬ 
trodes  concluded  from  Fig.  6  (about  120  s  and  370  s)  are  lower 
than  theoretical  values.  Increasing  lithium  ion  diffusivity  in  solid 
phase  may  decrease  the  polarization  and  improve  the  relaxation 
performance  of  battery. 

Besides,  Fig.  7  shows  the  potential  in  solid  phase  also  changes  at 
pulse  test.  The  potential  decreases  with  time  as  a  typical  saw-tooth 
behavior,  corresponding  to  the  increase  of  lithium  ion  concentra¬ 
tion  on  the  surface  of  positive  particle  (Fig.  6c  and  d).  Within  200  s, 
the  voltage  decrease  at  boundary  5  that  is  larger  than  that  at 
boundary  4  due  to  the  ohmic  overvoltage  contributed  from  the 
thickness  of  positive  electrode.  The  relaxation  of  the  potential  after 
current  interruption  occurs  in  two  stages.  When  the  current  is  cut 
down  at  200  s,  the  decrease  is  also  shut  down  immediately.  And  the 
potential  increases  by  0.04  V  quickly  with  a  timescale  less  than 
0.2  s,  which  is  caused  by  the  disappearance  of  ohmic  polarization. 
Duo  to  the  diffusion  of  lithium  ions,  the  potential  increases  slowly 
until  equilibrium  is  attained. 


4.2.  Dynamic  evolution  for  electrochemical  reaction 

In  order  to  observe  the  dynamic  evolution  in  the  electrodes,  the 
local  current  density,  surface  concentration  and  the  state  of  charge 
are  simulated.  The  curves  of  these  properties  at  different  time 
during  discharge  are  shown  in  Figs.  8  and  9.  In  the  negative  elec¬ 
trode  (Fig.  8a),  the  local  current  density  does  not  change  in  regu¬ 
larity.  For  the  positive  electrode,  the  local  current  density  at 
Boundary  4  increases  severely  at  the  beginning  of  discharge.  A  peak 
is  seen  at  1220  s  in  the  local  current  density  curve,  and  it  moves 
from  the  separator-positive  electrode  interface  to  the  positive 
electrode-current  collector  interface.  This  indicates  large  electro¬ 
chemical  reaction  rate  near  the  separator-negative  electrode 


Fig.  5.  Spatial  and  temporal  distribution  of  Li+  concentration  in  electrolyte. 
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Fig.  6.  Li+  distribution  along  particle  radius  during  pulse  tests  (boundary  2,  3,  4  and  5). 


interface.  This  result  is  different  from  the  report  of  Wang  and  his 
co-works  [22]  in  which  the  peak  moved  from  the  current  collector 
to  the  cathode-separator  interface  as  discharge  proceeded.  This  is 
because  the  peak  location  depends  upon  the  comparison  between 
the  effective  conductivity  of  the  solid  phase  (kff)  and  that  of  the 
liquid  (k|ff)  [51].  When  kff  <  /<|ff,  the  peak  appears  near  the  cur¬ 
rent  collector  and  moves  from  the  current  collector  to  the  cathode- 
separator  interface  as  discharge  proceeds.  On  the  contrary,  when 
kf  >  /<|ff,  the  peak  appears  near  the  cathode-separator  interface 
and  moves  from  the  cathode-separator  interface  to  the  current 
collector  as  discharge  proceeds.  In  our  work,  kff  is  larger  than 
therefore,  the  surface  concentration  near  the  separator-positive 
electrode  is  increased  rapidly  to  the  maximum  concentration  of 
22,806  mol  m-3  as  shown  in  Fig.  9a,  and  the  positive  particles 
approximately  reach  the  fully  state  of  charge  in  a  short  time  (about 
1800  s)  (Fig.  9b). 


The  state  of  charge  (SOC)  is  defined  by  the  quotient  of  the  sur¬ 
face  concentration  divided  by  the  maximum  concentration  of 
22,806  mol  m-3.  Therefore,  the  surface  concentration  and  the  state 
of  charge  distributions  at  different  time  represent  the  same  as 
shown  in  Fig.  9.  The  particles  near  to  the  separator-positive  elec¬ 
trode  interface  firstly  approach  the  fully  discharge  state.  As  a 
consequence,  the  utilization  of  active  materials  is  not  the  same.  As 
the  curve  shown  in  Fig.  8b,  the  peak  of  the  local  current  density  at 
the  end  of  discharge  has  not  reached  the  positive  electrode-current 
collector  boundary  at  t  =  3577  s.  Correspondingly,  the  surface 
concentration  in  the  region  close  to  the  current  collector  also  has 
not  reached  the  maximum  concentration.  In  addition,  the  local 
state  of  discharge  is  far  less  than  1,  which  indicates  that  the  active 
material  is  not  adequately  utilized.  That  is  because  the  stop  con¬ 
dition  is  determined  by  the  cut-off  voltage  in  both  simulation  and 
practice.  When  the  cut-off  voltage  is  reached,  the  utilization  of 
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and  5). 


active  materials  may  not  be  100%.  To  optimize  battery  design  and 
achieve  better  performance,  it  is  essential  to  choose  appropriate 
values  for  some  design  parameters. 

4.3.  Heat  generation  and  thermal  behavior 

Understanding  the  thermal  behavior  of  these  heat  generation 
sources  would  offer  much  valuable  information  to  develop  the 
battery  thermal  management  strategies.  There  are  several  heat 
generation  mechanisms,  as  shown  in  Fig.  10.  The  overall  heat 
generation  of  the  battery  discharge  is  consisted  of  electrochemical 
reaction  heat,  ohmic  heat  and  active  polarization  heat  produced  in 
the  negative  electrode,  electrolyte  and  positive  electrode.  The  heat 
generations  of  negative  and  positive  electrodes  are  complex,  as 
shown  in  Table  1.  The  ionic  ohmic  and  electrochemical  interface 
heats  are  dominant  contributor,  nearly  100%  of  the  total  heat 
generation.  The  ionic  ohmic  heat  from  negative  electrode  is  positive 
and  nearly  constant  through  the  whole  discharge  except  0-25  s, 
which  is  mainly  due  to  the  rapid  change  of  the  ionic  ohmic  at  early 
discharge.  An  S-shape  change  of  the  electrochemical  interface  heat 
is  observed  from  negative  to  positive  as  a  function  of  time.  Its  shape 
resembles  a  mirror  image  of  the  entropy  change  versus  SOC  curve 


Fig.  8.  The  local  current  density  distributions  at  different  time  during  the  discharge  at  1C  rate,  (a)  The  negative  electrode  side  and  (b)  the  positive  electrode  side. 


Fig.  9.  The  state  of  charge  distribution  at  different  time  during  the  discharge  at  1C  discharge. 
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of  the  LixC6  anode  in  Fig.  2a.  Due  to  the  electrochemical  interface 
heat,  the  total  heat  generation  rate  of  negative  electrode  also  ex¬ 
hibits  an  S  shape. 

As  shown  in  Fig.  10b,  the  value  of  ionic  ohmic  heat  from  positive 
electrode  is  also  positive  through  the  whole  discharge  progress.  But 
it  appears  a  large  range,  which  illustrates  that  the  ionic  ohmic  of 
positive  electrode  is  more  variable  than  the  negative  electrode. 
With  the  combined  effects  of  ionic  ohmic  heat  and  electrochemical 
interface  heat,  the  total  heat  of  the  positive  electrode  does  not 
exhibit  the  shape  as  the  entropy  of  the  positive  electrode  change 
versus  SOC  curve  of  the  LiFeP04  cathode  in  Fig.  2b.  Compared  to  the 
heat  generation  of  porous  electrodes,  the  heat  generation  of  the 
current  collectors  and  separator  are  tiny  and  can  be  neglected,  as 
shown  in  Fig.  10c  and  d. 

As  introduced  in  the  part  of  energy  balance,  the  heat  generation 
can  be  divided  into  two  types,  irreversible  heat  and  reversible  heat. 
The  heat  generation  of  the  current  collectors  and  separator  shown 
in  Fig.  10c  and  d  is  irreversible  joule  heat.  Fig.  11  plots  the  irre¬ 
versible  heat,  reversible  heat  and  the  total  heat  of  porous  electrodes 
and  the  battery  system.  Comparing  Eqs.  (19),  (24)  with  Eq.  (26),  we 
obtain  the  relation  as  Eq.  (45),  which  means  that  the  electro¬ 
chemical  interface  heat  of  porous  electrodes  comprises  reversible 
heat  generation  and  active  heat  generation.  The  curve  of  the 
reversible  heat  in  Fig.  11a  almost  has  the  same  shape  with  the  curve 
of  electrochemical  interface  heat  in  Fig.  10a,  except  the  time  range 
from  1450  s  to  discharge  end.  From  the  value  and  shape  difference, 


we  can  deduce  that  the  active  heat  generation  of  negative  electrode 
is  about  5  x  104  W  m-3  and  have  a  server  increase  from  1450  s  to 
the  discharge  end.  And  the  active  heat  generation  of  the  positive 
electrode  remains  a  stable  value  of  2  x  104  W  m-3.  The  active  heat 
generation,  also  called  polarization  heat  generation,  is  directly 
dependent  on  the  overpotential,  so  we  can  conclude  that  the 
variation  of  the  overpotential  in  negative  electrode  is  severer  than 
in  positive  electrode,  especially  at  the  last  stage  of  discharge.  The 
negative  electrode  plays  an  important  role  for  the  contribution  to 
the  overpotential  in  a  complete  battery.  The  total  heat  generation  of 
the  negative  electrode  is  slightly  higher  than  that  of  the  positive 
electrode.  While  the  thickness  of  the  positive  electrode  (70  pm)  is 
more  than  twice  that  of  the  negative  electrode  (34  pm),  so  the  heat 
generation  of  a  complete  battery  come  mostly  from  the  positive 
electrode,  which  is  the  reason  why  the  shape  of  the  total  heat 
generation  is  similar  to  that  of  the  positive  electrode,  as  shown  in 
Fig.  lib  and  c. 


Qi,(l)  —  ^a,iJloc,i  +  01, i  02, i 

(n  _  c  .  TdU{\ 

—  (  *2rea  —  *->a,i/loc,i  *  J 

+  [Qact  =  Sa,i/loc,  i  (01,i  —  02,  i  —  U‘i)] 


(45) 


As  shown  in  Fig.  11,  the  reversible  heat  accompanied  endo¬ 
thermic  process  and  exothermic  process,  and  the  irreversible 
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Fig.  10.  Time  history  of  heat  generation  by  various  mechanisms  during  discharge  at  2C  rate,  (a):  negative  electrode;  (b):  positive  electrode;  (c):  current  collectors;  (d):  separator. 
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Fig.  11.  Time  history  of  reversible  and  irreversible  heat  generation  during  discharge  at  2C  rate,  (a):  Negative  electrode;  (b):  positive  electrode;  (c):  complete  battery. 


heat  is  exothermic  progress  and  increases  with  increasing  the 
discharge  rate  [52].  Their  combination  determines  whether  the 
discharge  process  is  exothermic  or  not.  Fig.  12  shows  the  total 
heat  generation  of  complete  battery  at  different  discharge  rate. 
The  blue  (in  the  web  version)  area  stands  for  the  endothermic 
heat,  the  red  area  represents  the  exothermic  heat  with  an  equal 
amount  of  the  blue  area,  and  the  yellow  area  represents  the  net 
heat.  At  low  discharge  rate  like  0.5C,  the  discharge  electric  cur¬ 
rent  is  so  small  that  the  polarization  is  not  serious,  which  leads 
smaller  exothermic  amount  of  irreversible  heat  than  the  endo¬ 
thermic  amount  of  reversible  heat.  Hence,  there  are  three  colored 
areas  existing  in  Fig.  12a,  indicating  the  discharge  process  of  the 
LiFeP04  battery  is  accompanied  exothermic  and  endothermic 
processes  at  low  rates.  The  same  phenomena  still  exists  in  1C 
discharge  process.  But  when  the  rate  is  increased  to  2C  and 
above,  only  yellow  area  is  left,  indicating  only  exothermic  phe¬ 
nomena  existing.  This  phenomenon  is  in  a  good  agreement  with 
the  experimental  work  of  Liubin  Song  and  his  coworkers  [53,54]. 
Fig.  13a  shows  evolutions  of  the  average  battery  temperature 
during  1C  discharging.  The  battery  is  surrounded  by  natural  air  at 
25  °C.  At  the  initial  stage  of  discharge,  the  average  temperature 
increases  slowly.  But,  the  average  temperature  increases  to  about 
37  °C  at  the  end  of  discharge.  Fig.  13b  is  the  contour  plot  of  the 
battery  temperature  distribution  at  the  end  of  1C  discharging. 
Hot  area  is  distributed  at  the  axial  core  of  the  battery.  Due  to 
larger  thermal  conductivity  in  the  axial  direction,  the  tempera¬ 
ture  distribution  in  this  direction  is  uniform.  The  main  temper¬ 
ature  variation  is  in  the  radial  direction  with  a  value  about  1.5  °C. 
Basing  on  the  comparison  between  the  heat  generation  of  5C  and 
1C  in  Fig.  12  and  the  temperature  contour  plot  at  1C  in  Fig.  13,  the 
battery  temperature  will  reach  more  than  60  °C  at  5C  discharging 
rate.  Appropriate  rates  can  be  used  to  avoid  the  internal  heat 
accumulation  of  batteries. 

All  computations  are  carried  out  on  a  DELL  PRECISION  T1650 
Workstation  with  two  quad-core  processors  (Intel  Core  i7-3770, 
3.4  GHz,  with  a  total  of  eight  processor  cores)  and  a  total  of  16  GB 
random  access  memory  (RAM). 


5.  Conclusion 

In  this  paper,  an  electrochemical-thermal  model  based  dy¬ 
namic  materials  response  for  lithium  iron  phosphate  battery  is 
developed  by  employing  the  comprehensive  dynamic  parameters 
in  thermodynamics  and  kinetics.  The  current  collectors  are 
considered  in  the  model.  This  model  is  validated  in  aspects  of 
electrochemical  performance,  thermal  performance,  which  is  in  a 
good  agreement  between  the  simulated  results  and  experimental 
results. 

The  pulse  tests  show  that  the  self-adjusting  period  needed  for 
lithium  ion  concentration  in  electrolyte  to  attain  equilibrium 
after  relaxation  is  affected  by  the  discharge  rate.  The  greater 
discharge  rate  is,  the  longer  self-adjusting  period.  Because  the 
rate  of  lithium  ions  deintercalating  from  the  negative  porous 
electrode  is  larger  than  that  of  intercalating  into  the  positive 
electrode,  after  enough  relaxation,  the  lithium  ion  concentration 
in  electrolyte  shows  an  augmentation  compared  to  initial  level. 
And  the  self-adjusting  period  needed  for  lithium  ion  concen¬ 
tration  in  particle  to  attain  equilibrium  is  longer  than  that  for 
electrolyte.  The  existence  of  critical  makes  the  dynamic  evolu¬ 
tion  for  lithium  ion  concentration  in  solid  phase  much  compli¬ 
cated.  Moreover,  in  constant  current  discharge  processes,  the 
electrochemical  reaction  rate  is  location-related,  and  the 
maximum  rate  moves  from  the  separator-positive  electrode 
interface  to  the  positive  electrode-current  collector  interface, 
resulting  the  utilization  of  active  material  is  not  uniform.  At  the 
end  of  discharge,  there  is  still  a  part  of  the  active  material  has 
not  been  adequately  utilized.  The  discharge  process  of  LiFePCH 
battery  accompanied  exothermic  process  and  endothermic  pro¬ 
cess  at  a  low  rate.  With  the  rate  increasing,  the  endothermic 
phenomena  disappears  gradually,  the  heat  production  rate  and 
the  enthalpy  change  during  the  discharge  process  of  LiFePCH 
battery  increase.  When  the  rate  is  increased  to  2C,  only 
exothermic  phenomena  exists.  These  suggest  that  appropriate 
rates  for  batteries  should  be  chosen  in  order  to  avoid  safety 
problem. 
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Fig.  12.  The  total  heat  generation  of  complete  battery  at  different  discharge  rate. 


Fig.  13.  (a)  Evolutions  of  the  average  battery  temperature  during  1C  discharging  and  (b)  the  contour  plot  of  the  battery  temperature  distribution  at  the  end  of  1C  discharging. 
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